THz emission from a stacked coherent flux-flow oscillator: non-local radiative 
boundary conditions and the role of geometrical resonances. 



O 
(N 



l> 



V. M. Krasnov 
Department of Physies, Stockholm University, AlbaNova University Center, SE- 10691 Stockholm, Sweden 

(Dated: May 18, 2010) 

I derive simple non-local dynamic boundary conditions, suitable for modelling of radiation emis- 
sion from stacked Josephson junctions, and employ them for analysis of flux-flow emission from 
intrinsic Josephson junctions in high-Tc superconductors. It is shown that due to the lack of Lorenz 
contraction of fluxons in stacked junctions, high quality geometrical resonances are prerequisite for 
high power emission from the stack. This leads to a dual role of the radiative impedance: on the 
one hand, small impedance increases the efflciency of emission from the stack, on the other hand, 
enhanced radiative losses reduce the quality factor of geometrical resonances, which may decrease 
the total emission power. Therefore, the optimal conditions for the coherent flux-flow oscillator are 
achieved when radiative losses are comparable to resistive losses inside the stack. 

PACS numbers: 74.72.Hs, 74.78.Fk, 74.50.+r, 85.25.Cp 
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Creation of a compact, high power THz source remains 
a difficult technological challenge, colloquially known 
as "the THz gap" [l]. Flux-flow oscillators (FFO's), 
based on regular motion of quantized vortices (fluxons) 
in Josephson junctions, can generate tunable THz radi- 
ation with a remarkable linewidth 10~^^, albeit with a 
low emission power < IfiW ^. The power could be 
greatly enhanced by coherent phase-locking of several 
coupled FFO's @, |j]. The required coupling is strongest 
in stacked, atomic scale intrinsic Josephson junctions 
(IJ Js) , naturally formed in layered high temperature su- 
perconductor Bi2Sr2CaCu208-i-a; (Bi-2212). IJJs allow 
simple integration of a large number of almost identical 
stacked Josephson junctions. Furthermore, a large super- 
conducting energy gap in Bi-2212 jSi] facilitates operation 
in the important THz frequency range (up to ^ lOTHz). 
Therefore, IJJs are actively studied both theoretically 6- 
[ll| and experimentally |12| - [l7j as possible candidates for 
realization of a coherent THz oscillator. 

Proper radiative boundary conditions are essential for 
analysis of the coherent FFO. For a single junction this 
can be done by introducing an appropriate radiation 
impedance Z [l8| , connecting local ac-components of elec- 
tric and magnetic fields at the junction edges. 



Z = Z{u) = Eac/Ha 



(1) 



The emission power is Prad = Vac/^Zi where Vac — toEac 
is the ac-voltage, to is the junction barrier thickness and 



Rz = ito/w)Z 



(2) 



is the radiative resistance of a single junction, w is the 
width of the junction. However, the stacked FFO can not 
be described by a single Z because the effective radiative 
resistance depends not only on the geometry of the stack, 
but also in a crucial way on the collective fluxon config- 
uration. For example, motion of the triangular fluxon 
lattice corresponds to Z ^^ oo because it results in out- 
of-phase oscillations in neighbor junctions Ei — — i?i+i, 
leading to destructive interference and negligible emis- 
sion. On the other hand, motion of the rectangular lattice 



corresponds to small Z because it results in in-phase os- 
cillation Ei = Ei+i, leading to constructive interference 
and coherent enhancement of the emission power oc N'^, 
where N is the number of junctions in the stack. Such 
behavior is caused by the essentially non-local nature of 
magnetic induction [19], which is the consequence of the 
inductive coupling of junctions in the stack. Therefore, 
the relation between local Eac and non-local Hac can not 
be described by Eq. ([T]). 

In this letter I derive simple non-local dynamic radia- 
tive boundary conditions for stacked Josephson junctions 
and employ them for analysis of flux-flow emission from 
IJJ-based stacked FFO. It is shown that due to the lack 
of Lorenz contraction of fluxons in stacked junctions 'l£|, 
high quality Q 3> 1 geometrical resonances are prerequi- 
site for high power emission from the stack. This leads 
to a dual role of the radiative impedance: on the one 
hand, small Z increases the emission efficiency, but on 
the other hand, enhanced radiative losses reduce Q of 
geometrical resonances, which leads to reduction of the 
emission power. The maximum power is achieved when 
radiative losses are equal to internal losses in the stack. 

Neumann (static) boundary conditions are most widely 
used for numerical modelling of Josephson junctions: 






(3) 



Here ip is the Josephson phase difference, x — 0,L are co- 
ordinates of junction edges, B is the applied dc-magnetic 
field, Xj is the Josephson penetration depth and A* is 
the effective magnetic thickness of the junction [19[. As 
emphasized in Ref. |18| , Neumann boundary conditions 
are non-radiative because they assume HadO, L) — 0. 

Radiative (dynamic) boundary conditions should ac- 
count for the finite Hac- For a single junction they can 
be easily written with the help of Eq. p|) [TTI,[l8i: 
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(4) 



Here plus and minus signs correspond to x = and L, 
respectively, because the direction of emission is opposite 
at opposite edges of the junction. 

However, Eq.([3]) is not directly applicable for stacked 
junctions due to the principle difference in electrodynam- 
ics, which is local in single, and non-local in stacked junc- 
tions. The electric field Eac is local in each junction be- 
cause the Debye screening length is always smaller than 
the electrode thickness d, even for IJJs. However, the 
magnetic field in inductively coupled stacked junctions 
is non-local and is created and shared cooperatively by 
the whole stack [1^ . The non-locality of Hac is particu- 
larly dramatic for IJJs, for which the London penetration 
depth Xab — 200 nm is much larger than d < 1 nm. 

To derive the proper non-local radiative boundary con- 
ditions for stacked Josephson junctions, lets note that the 
local relation Eq.(IT]) is valid outside the stack, however, 
here the electric field is the result of interference of elec- 
tric fields from all junctions: 



JV 



Eac — 2_^ ^i 



(5) 



The net emission power from one edge of the stack is 

N 
Prad = wtoHac ^ ^^ = wtoSU Z . (6) 

i=l 

Radiative losses are associated with additional currents 
flowing through edges of the stack: 



A/rad = wEac/Z 



(7) 



Those displacement-like currents should be added in the 
numerical scheme for the sake of energy conservation. 
Equations (I4I5I7I) together form the final non-local dy- 
namic boundary conditions for stacked Josephson junc- 
tions. For the out-of-phase state, Ei = — i?i+i, Eac = 0, 
they reduce to the non-radiative Neumann condition 
Eq.®. For the in-phase state, E^ — i?j+i, Eac = NE^, 
they lead to coherent power amplification ex A^^, Eq.®. 
To calculate flux-flow emission, the non-local bound- 
ary conditions were implemented into the coupled sine- 
Gordon equation [J|, which describes electrodynamics 
of inductively coupled stacked Josephson junctions (see 
Refs.|19l. l2q | for details of the formalism and the Sup- 
plementary [2ll | for technical details). Simulations were 
made for TV = 10 identical IJJs with L = 4Aj ~ 2.8 /im, 
w = 1.4 /im, the stacking periodicity s = 1.55 nm, the 
electrode thickness d = 0.6 nm, the ratio of the barrier 
thickness to the relative dielectric constant io/cr — 0.11 
nm, Xab — 200 nm, and the critical current density Jcq = 
1050 A/cm^. Those parameters correspond to Xj ~ 
0.7/im, the Josephson plasma frequency top ~ 6.4 x 10^^ 
1/s, and the slowest velocity of light cat ~ 3.2 x 10^ m/s, 
typical for small Bi-2212 mesa structures [l^jlS^l- The 
damping parameter a = {upRC)^^ — 0.05, corresponds 
to junction resistance R ~ 100 fl and the c— axis resis- 
tivity pc ~ 25 r^cm, equal to large bias tunnel resistivity 
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FIG. 1. (color online). Simulated a.) I — V characteristics 
and b) normalized radiative powers for negligible Z — 10"^^ Q 
(dark blue) and moderate Z — 2 x 10* fl (red) emission. A 
variety of Fiske steps is seen. Mode numbers {m, n) for the 
most prominent steps are indicated. Note that non-radiating 
even-n steps (dashed vertical lines) are not affected by Z, 
while radiating odd-n steps (solid lines) vanish for small Z. 



of IJJs 23]. However, the quasiparticle (QP) resistivity 
at small bias and low T can be up to two orders of mag- 
nitude larger ^]. Therefore, a = 0.05 represents the 
highest limit of QP damping in IJJs. The correct a is 
important for numerical modelling because it determines 
Q oc 1/a in the absence of radiative losses. Simulations 
were also performed for smaller a down to 0.005, which 
show qualitatively similar behavior [17| , but require sub- 
stantially longer integration times. However, for larger 
a = 0.1 the number and amplitudes of accessible geo- 
metrical resonances were significantly reduced. 

Fig. 1 a) shows simulated I ~ V characteristics for 
two impedance values, corresponding to negligible Z = 
10^2 Q (dark blue) and moderate Z = 2 x W n (red 
curve) radiation losses, at ^/^o — 3, where $ = BLs is 
the flux per junction. Data for Z = 10'^^ fi are similar to 
simulations with non-radiative Neumann boundary con- 
ditions [a, 0, Il3 • A large variety of Fiske steps is seen. 

Geometrical (Fiske) resonances are caused by resonant 
excitation of phonon-like collective fluxon lattice vibra- 



tions T?l in the flux-flow state. In stacked junctions such 
fluxon phonons are two-dimensional and are character- 
ized by two wave numbers km — inn/L {m — 1, 2, 3, ...) 
in plane, and qn — im/Ns [n — 1,2, ...N) in the c-axis 
direction [25]. Unlike Josephson plasma waves, fluxon 
phonons have a linear dispersion relation at low frequen- 
cies. Fluxon phonons with the wave number (/„ propagate 
with the in-plane velocity c„ [25|. The slowest, cat, cor- 
responds to the out-of-phase, qn = i^/s, and the fastest, 
Ci, to the in-phase, gi = tt/Ns, mode. Fiske steps oc- 
cur when the ac- Josephson frequency coincides with the 
frequency of one of the cavity modes (m, n) |l7[: 



Vm 



<^omCn/2L, (m = 1,2,3. 



1,2...N). (8) 



In Fig. 1, voltages of the most prominent phase- locked 
Fiske steps NVm,n are marked by vertical lines together 
with mode numbers (m,n). For the non-radiative case 
Z = 10"^^ ft, Fiske steps with both even and odd n are 
equally represented. However, for Z = 2 x 10* U we 
observe a clear difference between even and odd n reso- 
nances: even-n steps, e.g. (1,10), (3,10), (5,10), (6,4), 
(7,4), (6,2) and (7,2) are practically the same as for 
Z = 10"^^ ft, but all odd-n steps are strongly reduced. 
For example, steps (5,3), (5,1) and (6,1) are well devel- 
oped for Z = 10^2 ^^ bu^ noi- present for Z = 2 x 10^ n. 

The observed difference between even and odd-n steps 
for Z — 2 X 10^ il is due to appearance of significant 
radiative losses. For the considered case of even N = 10, 
all even-n resonances are interfering destructively and 
therefore are non-emitting. To the contrary all odd-n 
resonances are emitting progressively stronger with de- 
creasing the mode number n. This is most clearly seen 
from the voltage dependence of the normalized emission 
power, shown in Fig. lb). It can be seen that for 
Z = 10'^^ ri even-n Fiske steps (dashed lines) correspond 
to minima, while odd-n steps (solid lines) to maxima 
of emission. As expected, the largest (by two orders of 
magnitude) emission occurs for in-phase resonances (5,1), 
(6,1) with the largest Fiske step amplitudes, see Fig. 1 
a). They correspond to the velocity matching condition 
m ~ 2int[$/$o], as discussed in Ref. 17]. The Fiske 
step (5,3) has larger amplitude but smaller emission be- 
cause only one third of the stack is emitting coherently 
while the rest of the junctions are interfering destruc- 
tively. Also odd-n Fiske steps (1,3), (3,3), (6,3), (7,3), 
(4,1), (7,1) with small amplitudes, see Fig. 1 a), produce 
clear emission peaks for Z = 10^^ fi. 

Remarkably, all these emission peaks are absent for 
Z = 2 X W^ n. Furthermore, the strongest emission 
peaks (5,1), (6,1) at Z = 10^^ il are replaced by minima 
for Z = 2 X 10** il. Apparently, radiative losses have com- 
pletely damped in-phase resonances for Z = 2 x 10^ SI. 
This clearly illustrates that the quality factor is playing 
a central role in operation of the coherent FFO. 

The role of the output impedance on the FFO power 
can be understood from the following simple model. The 
oscillating amplitude at the resonance is Ei = EqQ = 
EoujReffC, where Eq is the amplitude out of the reso- 
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FIG. 2. (color online). Simulated emission characteristics 
(symbols) as a function of the radiation impedance for three 
voltages, marked in Fig. 1. A: at the out-of-phase (1,10) reso- 
nance, B: in the incoherent state, and C: at the in-phase (6,1) 
resonance. Lines represent analytic expressions Eqs. (|10llip 
with Rrad indicated in b). It is seen that the emission power 
is decreasing both at small and large Z and that the max- 
imum power is achieved at Rrad = R, indicated by arrows 
in a). Note that the same stack has different values of Rrad, 
depending on collective fluxon dynamics in the stack. 



nance and Reff is the effective damping resistance due 
to both internal QP (i?) and radiative {Rrad) losses: 



^e// - [R^ 



R 



'rad\ 



(9) 



The total emission power from one side of the junction is 

R Rrad 



Prncl — EniO tnW , ,^ . 

° (R + RradY 

It is useful to consider the quantity 

PradZ = 7(l+i?/Kad)"', 



(10) 



(11) 



where 7 is constant for a given resonance uJm,n- 

As already discussed above, see Eqs. (|5I6I) . radiative 
losses in stacked Josephson junctions depend essentially 
on the collective, rather than local, fluxon dynamics. 
They can be described by the "coherence factor" : 



Cp = \Eac/Ei\ , 



(12) 



which reflects the efficiency of emission in comparison to 
a single jmiction with the same electric field amplitude. 
The effective radiative resistance (per junction) is then: 



Rrad — Rz/Cp- 



(13) 



Fig. 2 shows a) the emission power (from one edge) 
and b) the product PradZ, as a function of Z for three 
voltage levels, marked by arrows in Fig. 1: A - the co- 
herent, phase- locked state at the lowest out-of-phase ge- 
ometrical resonance 14, lo — 1-1 mV; B - an incoherent, 
chaotic state at F ~ 9 mV at which no clear resonances 
are observed and individual junctions are having different 
voltages; C - the coherent in-phase resonance Ve,! — 47 
mV. Symbols and lines represent numerical simulations 
and analytic expressions, Eqs. (jlOllip . respectively. It is 
seen that at large Z, PradZ is constant. In this case radia- 
tive losses are small compared to QP losses, Rrad ^ R^ 
the amplitude Eac is independent of Z and Prad is in- 
creasing oc 1/Z with decreasing Z, see Eq.®. However, 
for Rrad ^ Ri ths quality factor starts to decrease, re- 
sulting in reduction of Eac and the emission power. 

According to Eq. (|T0)) . the maximum of Prad is 
achieved when Rrad — R- However, as follows from Eq. 
([T^. Rrad for stacked junctions is not constant, but de- 
pends on the collective state of the stack, described by the 
coherence factor Cf- For the in-phase state C, Cf = N 
and Rrad = Rz /N . The corresponding Prad and PradZ 
calculated from Eqs. (llOllip are shown by solid lines 
in Figs. 2 a) and b) and describe well the numerical 
data. In this case the maximum of Prad is achieved at 
large Z — RNw/to, marked by the blue arrow in Fig. 2 
a). The incoherent /chaotic state B, is well described by 



Eqs. (llOllip with iV— times larger Rrad — Rz (dashed 
lines), because junctions in the stack act as individual 
single junctions, Cf — 1. Finally, as expected, the out- 
of-phase state corresponds to the largest Rrad which is, 
however not infinite, but fairly well described by yet an- 
other A^— times larger Rrad — NRz, Cf — 1/A^, repre- 
sented by the dashed-dotted lines. 

To conclude, strong emission from the stacked coherent 
fiux-fiow oscillator requires both large coherence factor 
Cf and large oscillating amplitude Eac- For a single junc- 
tion FFO, Eac can be greatly enhanced due to Lorentz 
contraction of the fiuxon [26| at the velocity matching 
condition [2]. However, such mechanism is absent in 
stacked Josephson junctions due to the lack of Lorentz 
invariance in the system [191 . Therefore, large Eac can be 
achieved only in the presence of a high quality resonance. 
Importantly, high quality geometrical resonances can im- 
pose their order on the fiuxon lattice. E.g., fast in-phase 
fiuxon phonons may stabilize the rectangular fiuxon lat- 
tice with large Cf = A^ even at low fiuxon velocities 
[ITj . In addition, the resonance reduces the linewidth of 
emission Sf/f oc l/Q- Therefore, high speed geometri- 
cal resonances are particularly important for realization 
of the coherent stacked FFO. As shown above, optimal 
operation of the coherent FFO is achieved when radia- 
tive losses are equal to internal resistive losses. Finally, 
it was shown that radiative losses lead to the asymme- 
try between even/odd-n Fiske steps in / — T^s of stacked 
Josephson junctions. This is consistent with resent ex- 
perimental observations for Bi-2212 mesas [13] and can 
be viewed as indirect evidence for significant coherent 
fiux-fiow emission from intrinsic Josephson junctions. 
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